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In this paper we compare two loop Chiral Perturbation Theory (ChPT) calculation of pion-pion 
scattering with the unitarity second order correction to the current algebra soft-pion theorem. It is 
shown that both methods lead to the same analytic structure for the scattering amplitude. Current 
algebra; Unitarity. 
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I. INTRODUCTION 



In the early sixties, many results for low energy meson physics have been derived from the assumptions of a 
local chiral SU(2) xSU(2) algebra of vector and axial vector current densities together with the partial conservation 
hypothesis (PCAC) relating the derivative of the axial vector current to the pion field. By itself, the PCAC relation 
can not help one to extend the applicability of current algebra method to the energy corresponding to meson-meson 
resonances. Nevertheless, a method was invented to obtain results for process were mesons are not "soft particles" . 
f*"*- \ We are referring to the hard meson methods of current algebra 0. Even ignoring the underlying theory, the chiral 
current algebra implies a set of Ward identities and the method consists in solving the system of Ward identities under 
the general principles of analyticity, crossing and elastic unitarity. 

In 1979, Weinberg suggested that it is possible to summarize these previous results in a phenomenological Lagrangian 
which incorporates all the constraints coming from chiral symmetry of the underlying theory 0]. In a set of very 
important and fundamental papers, Gasser and Leutwyller have developed Chiral Perturbation Theory (ChPT) which 
O i| allows one to compute many different Green functions involving low energy pions || . One of the main obstacles for the 
applications of ChPT to high energies lies in the issue of unitarity. Several methods have been proposed to extend the 
range of energies where ChPT could be applied. These methods include the Pade expansion 0, the inverse amplitude 
method || and the introduction of fields describing resonances . 

Let us remember, however that, in order to treat pion-pion scattering amplitude obtained by the hard meson method 
based on the Ward identity technique, one of us have introduced the constraints of elastic unitarity for partial waves 
0. We will call this approach by unitarization program of current algebra (UPCA). It was applied to obtain first 
order corrections (QUI) to the soft-pion 7T7t Weinberg amplitude || as well as to calculate second order corrections 
(QU2) to the referred amplitude M. 

The ChPT calculation of pion-pion scattering to one loop was performed in Ref. [Bj and, only recently, two loops 
calculation appeared in the literature [(Tof . Our aim is to compare UPCA and ChPT calculations. In a previous paper, 
we have shown that one loop ChPT is equivalent to QUI and in the present paper, we will compare QU2 amplitude 
with ChPT Lagrangian two loops calculation performed in Ref. Jl0|]. We will conclude that, as conjectured by one of 
us ||, the two approaches are equivalent. In the next section we will recall the comparison at one loop level and in 
the section three this comparison is extended to the two loops case. In the section four we present the conclusions. 
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II. FIRST ORDER CORRECTED BY UPCA AMPLITUDE AND ONE LOOP CHPT 



Let us remember the main points of the UPCA. The starting point in our derivation was an exact hard-pion 
expression for the correlation function of four currents, with the quantum numbers of the pion, in terms of three- and 
two-point functions. From this expression, by using vertex and propagator estimates, we could reobtain the so-called 
soft-pion Weinberg amplitude, namely: 



A CA (s,t,u) = ^(s-m 2 ). 



(1) 



The remaining of the amplitude reflects the difference between soft- and hard-meson result. In the UPCA one estimates 
the behavior of form factors and propagators at low energies and assumes that, for instance, for small values of its 
argument, the scalar pion form factor, Fx, is of the same order of magnitude as current algebra amplitudes near 
threshold. This can be obtained by setting, for x — m , Fi(x) ~ 1 + ff\x) + 0(e 2 ), for 7 = and 2, where 
the superscript "(1)" denotes the order e = m 2 /M 2 , M being of order of magnitude of vector meson masses. 

To construct unitarized amplitude, we have observed that current algebra gives real partial waves. The unitarization 
method must provide an imaginary part to the corrected partial wave. Thus, at the first order of the calculation, by 
the optical theorem, one must have 



Im $(s) = -L*(s) tf A (s) 2 
dZn 



where tf A is the soft-pion i partial wave, isospin /, Weinberg amplitude obtained from Eq. (fil) and 



Am 2 



a(s)- 

In the program, we work with exact total amplitude expression which follows from the Ward identities, and we 
use the implications of elastic unitarity for form factors and propagators in a peculiar way. For instance, from the 
relations valid for the scalar form factor and scalar propagator, namely: 



Im Fi{s) = ^- a(s) T/(s) Fx(s) and Im A/(s) = -J- a(s) |^/(s)| 2 , / = and 2, 
oZir oZn 

we obtain within the first order of the approximation: 



(2) 



-±-<t(s) tg A (s) and Im ^(.s) = ^-a(s), 

6ZTT oZit 



where t^ A is the current algebra isospin zero S-wave 7T7t amplitude, 



F 2 



(2s -m 2 ). 



Considering the known imaginary part of each function entering into the amplitude, the method consists in obtaining 
their real parts by the dispersion relation technique. To converge, dispersion integrals need subtraction which are 
model free parameters. They can be fixed by fitting experimental data. 

In this way, the first order corrected amplitude, QUI, derived in the context of UPCA [Q, can be written in the 
following form: 

F^l( S ,t,u) = i(2« - m 2 )^ 1 ^) - I (2m 2 - (s) + ^(s - 2M 2 ) 2 + 

i(2m 2 - t)^\t) + (s - u)*P(t) - - 2M 2 ) 2 + (*«-> u) 



with: 



^\x) = (2x-m 2 )(g(x)+a ), $£\x) = (2m 2 - x)(g(x) + a 2 ), 
$i } (x) - ^(x-4m 2 )g(x) - - (2x0, - 4m 2 g{0)) , 
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where: 

32 7r 2 5 ( S ) = (s-V) r dx- = a(s) In 44^ - ( 3 ) 

J 4m2 (x - 4m 2 ) (x-s) a(s) + 1 

On the other hand, the one loop elastic pion scattering obtained from ChPT Lagrangian || is: 

F* AW pT (s,t,u) = I( s 2 _m 4 )J( s )+ 

- [(t(t -u)- 2m 2 t + Am 2 u - 2m 4 ) J(t) + (*<-> u)] + 
[2(4 - 4/3)(s - 2m 2 ) 2 + (4 - 5/6)(s 2 + (t - uf)+ 

12m 2 s(£ 4 - 1) - 3(£ 3 + 44 - 5)m 4 ] /96tt 2 . (4) 

We can identify the function J(x) with 2 [ g(x) — g(0)) and we have verified that the polynomial coefficients of these 
functions are the same. We have then concluded that the above amplitude has the same analytical structure than 
QUI. Each approach have its free parameters: the model free parameters of QUI are £'s and a's linear combinations 
and the free parameters of ChPT are £'s linear combinations. In UPCA the free parameters are subtraction constants, 
inherent to dispersion relation technique, and in ChPT they come from tadpole graphs and tree graphs of order 0(p 4 ). 

From this comparison we have shown, in a recent letter, that one loop ChPT amplitude can fit experimental S- and 
P- waves up to the resonance region by adjusting only two parameters, namely 4 and 4 p| . 

In the next section, we will compare the second order corrected by UPCA amplitude, QU2, with two loops calculation 
|jTof . To do this, we will need one loop partial wave corresponding to the ChPT amplitude given above. For this, one 
expands the combinations with definite isospin in the s-channcl into partial waves: 

T I {s,t)=2,2-K ££L (2£ + 1)-Fi(cos0) t$(a), I = 0,1 and 2. 

Using Eq. (^), the resulting one loop P-wave amplitude is: 

= 7 ArAs~ 4m 2 ) 2 J(s) + l^( S 4- - m^-^W 



18F 4V ' w (_ 8 V2 16 J (s-4m 2 ) 

s 3 1 9 9 1 A m 6 \ L(s) 
— s wr + -sm 4 ' 



288 18 4 8 / er(s - 4m 2 ) 

- — . 1 ; 9N (s 3 + 37s 2 m 2 - 149.sm 4 + 120m 6 ) 
864 (s - 4m 2 ) V / 

c — zl 7 7 7 " 

' [(2£ 2 -24-l)s + 8m 



288 LV ' J I FS 2 

the isospin 7 = one loop S-wave is: 



4 1 o ) ( S ) = ^(2 S -- 2 ) 2 J>) + ^{^( i 



6 / (s — 4m 2 ) 



/ 7 2 5 2 25 4 \L s / 11 - 7 - 7 , 

s 2 s m 2 H m 4 — — + 4 H 4 s 

V144 18 48 / er V 144 72 96 1 



2 



n U 481\ 2 m 4 /ll- - 45- 9- 355 

— (4 + 4) — 4 sm 2 H — £1 + 74 4 — 4 H 

18 v y 4 432 / 18 V 2 16 4 16 



and the isospin 1 — 2 one loop S-wave is: 



2F 4V ; v ; F 4 tt 2 L 16 V 3 / (s-4m 2 ) 



ils 2 


1 




sm 


288 


9 








- —4 4 
18 








>- 



48 

5 \ 2 /I, 7 - 1 - 527\ 2 

s - — 4 H 4 + -4 sm 2 

192 J \36 36 8 864/ 

+ ^ + ^-^ + ^-i)m*}. (5) 

In these expressions we have included the contributions from J(4m 2 ) which is lacking in the expressions of partial 
waves given in the Sec. 2 of the Ref. lilt] and we have corrected an overall sign. 
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III. SECOND ORDER CORRECTED BY UPCA AMPLITUDE AND TWO LOOPS CHPT 



We have shown that the first order correction to the soft pion amplitude (QUI) is equivalent to one loop ChPT 
scattering amplitude and, in addition, we have given the tools for constructing the next order unitarity corrections 
(QU2) f§. The formula (3.10) of Ref. @ can be written as: 



F2AW 2 (s,t,u) = ±(2s-m 2 )^\s), 



±(2m 2 - s)<P {2 \s) 



^(2m 2 - t)^\t) + (s - u)^[ 2 \t) + {t^u) 



We can relate the above expression with the formula (3) obtained as a consequence of the Goldstone nature of the 
pion namely: 

F^Wj(s) = -Lt^( s )$f( s ), for 7= 0, 2 and F* W\(s) = j±-$?\s). 

OZTT 487T 

The functions Wi(s) are analytic except for a cut singularity at s > Am 2 . Their discontinuities are directly related 
with the discontinuities of the functions $/(s). 

We would like to emphasize that, the general structure of the UPCA solution comes from the Ward identity method 
and the hard meson technique implies that the amplitude is written in terms of form factors and propagators. We 
stress that the UPCA is based on the implications of elastic unitarity relations for form factors and propagators, and 
not for partial waves themselves. 

The consequences, of using Eq. (||), for instance for scalar form factors and propagators, to second order of the 
approximation, are: 

Im/f (s) = ^a(s) (EBt^(s)+t CA (s)BefP(s)), ImS^(s) = J-2Re/j 1) ( S ) 
with 1 = and 2. The vector form factor and vector propagator are obtained in a similar way. The functions 

(2) 

$) (s), constructed from form factors and propagators are then discontinued on the right hand cut as follows: 

Im$j 2) (x) = t^-ct(x) 2Re4 1} (a;), for I = 0, 1 and 2 

327T 

and Re t^ stands for the real parts of the functions in Eq. (Q). Using dispersion relation technique we obtain: 

$( 2) ( s ) = p I (s)Z(s) + q I (s)G(s)+n(s)g(s)+P I (s) for I = 0,1 and 2. (6) 

The polynomials pi(s), qi(s), rj(s), and Pi{s), for each value of total isospin I, are given in the Appendix. The 
function g(s) is given in (3)and 

(32 = ('-^ 2 )£-fc (l _ 4 : ( 4-,) (^)+"> 2 = lL(s)(L\s) + *>) 

where 

o-(s) - 1 



L(s) = In 



The amplitude is then written in terms of powers of L(s) and L(t) and contains the free parameters t\, 

^2 and the values of the subtraction constants ($/ and its derivatives at s — Am 2 ). 
On the other hand, ChPT amplitude calculated at two loops level Q with = 1 is: 

F*A% pT (s,t,u) = F( 2 \s)+G^(s,t) + G^(s,u), 
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with: 



FW(s) = J(s) 



1 



503 



16tt 2 V 108 



929 
"54 ' 



887 

~ 7 n i 



140 
~9~ 



+ 6i(4s-3) + 6 2 (s 2 +4s-4) 



--& 3 (8s 3 - 21s 2 + 48s - 32) + -& 4 (16s 3 - 71s 2 + 112s - 48) 



20s d 



19s 2 + 210s - 135 - — tt 2 (s - 4) 
16 



+ ^K 2 (s) (sir 2 - 24) + ±K 3 (a) (3s 2 - 17s + 9) , 



G< 2 )(s,<)= J(t) 



1 



16tH 



412 s , 9 

(t 

27 54 v 



+ -& 4 (2s(3t - 4)(t - 4) - 32i + 40r - lit ) 
6 



1 

36 
1 



Xi(i) 



174 + 8s - 10i 3 

1 + 4s + -^ 2 i(3s - 8) 
64 



72i 2 - 185* 



V216 


727 
108 


-4) 2 (3^ 
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-2s- 




4) (3s 



1571 \ 
108 ) 



9 

+-K 3 {t) (l + 3st - s + 3i 2 - 9t) + -i^ 4 (t) (4 



*)• 



In this expression J(s) = 2 (g(s) — g(0)), and: 
L 2 (s 



Ky = 

K 4 = 



(Wir 2 ) 2 
1 



, (16^ 2 ) 2 K 2 = a 2 L 2 (s) - 4, (16tt 2 ) 2 K 3 = — L 3 (s) + —L(s) - 

sa s a 1 



- K x + -K 3 

2 1 3 3 



1 - ir 2 
J + s ■ 



6 



16tt 2 



192vr 2 



(7a) 



(7b) 



Our strategy to compare AS~^^ with A^ pT was to expand them in terms of L(s) and L(t) and then to confront 
its coefficients. We have checked that they are the same. With respect to the polynomials we also realize that the 
structure are the same, but clearly the coefficients are differents because they have different origins. These polynomials 
include the model free parameters to be used in order to fit the available experimental data O] . 



IV. CONCLUSIONS 



Our aim is to compare Chiral Perturbation Theory (ChPT) calculations with the Unitarization Program of Current 
Algebra (UPC A). In previous works we have compared one loop ChPT with the first order corrected by UPCA 
results for pion-pion || and for kaon-pion scattering [|l3| and we have concluded that they lead to the same analytical 
structure for the amplitudes. 

The two loops calculation of pion-pion scattering only recently appeared in the literature [^0|. However, the tools 
for constructing second order corrections for the soft-pion current algebra result has been presented more than ten 
years ago Q. In the present paper we compare these results and we show that, as it was conjectured ||, they have the 
same analytical structure. We will shortly present the tools for constructing next order UPCA correction to kaon-pion. 
ChPT to two loops for Kn are not yet available, but we expected that it will be equivalent to the UPCA result too. 

In fact, the equivalence between the two approaches were expected. In the hard-meson method one starts from 
the chiral symmetric Ward identity exact result for the correlation function of four currents carrying pion quantum 
numbers. On the other hand, ChPT describes the low energy dynamics of fields realizing non-linearly chiral symmetry. 
Our unitarization program is based in the principles of analyticity, crossing and elastic unitarity which in turn are 
inherent to a field theory such as ChPT. 

In the framework of the Generalized Chiral Perturbation Theory, it was shown how to implement elastic unitarity 
strating from the one loop partial waves. That procedure leads to equivalent 0(p 6 ) ChPT amplitude ]1^] . However, 
the main difference from our unitarization procedure is that it uses the consequences of elastic unitarity for the form 
factors and propagators rather than for the amplitudes themselves. 
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Despite the fact that ChPT is a well stablished low energy effective theory for meson processes, we have shown here 
that the UPCA is a suitable alternative. One obstacle related to these two approaches is how to fix the free parameters 
which, in principle, are related to the parameters of the fundamental theory. At present, since this relationship is still 
lacking, the two approaches are both in same footing. In ChPT context we have shown that the one loop parameters 
4 and 4 can be fixed by fitting S- and P-wave tttt phase-shifts p~l| |. The final expression of pion-pion ChPT 
amplitude up to two loops diagrams has six parameters. However, the D-wave amplitude from pick up an 

imaginary part for s > 4 m 2 and we claim that the global fit of S-, P- and D-wave phase shifts will allow one to fix 
the new free parameters pa] . 



APPENDIX: 



The polynomials multiplying the functions g(s), G(s) and Z(s) in Eq. (Q) of second order correct by UPCA 
amplitudes are: 

2 , s f u i 7 - 17\ 2 (h- 5- 1- 95 \ 
7T 2 p (s) = —h + — 4 + -7T )s 2 - -4 + -4 - -4 - — s + 



72 36 48/ \9 9 2 72, 

11 - 7- 45- 9- 373 

— 4 4 -4 4 4 H , 

18 9 8 2 144' 

7T 2 pi(«) = 4^27^4 {( 6£ "i " ^"2 " 2 ) s 3 - (484 - 484 - 61) s 2 + 

(96^ - 964 - 197) s + 120} , 

^2 ( s) = ( }_- t + li + s 2 _ ( + l_j + \~ t + M_\ 
1 36 9 96/ I 18 18 4 4 144 i 



-4 4 -4 - -4 + -4 + — ■ 

9 9 8 218' 

q (s) = - 2 (50s 3 - 260s 2 + 303s - 36) , q 1 (s) = - 8 (6s 2 - 55s + 64) 



9(s-4) v /> 9( s _ 4 )5 

©(*) = - 9(g ^ 4) (10s 3 - 52s 2 + 93s - 72) , r (s) = 3(a ^ 4) (6s - 25) . 

ri(s) = (3s2 ~ 13s " 6) ' r2[s) = ~W^) (3s + 1} ■ 

The polynomial part of QU2, Pj, in Eq. (0) are written as: 

Pj(s) = A lS 2 + Bis + Ci, for I = and 2 
(s - 4) Pi = A lS 3 + B lS 2 + ds + D x 



with: 



= — £ + — * 2 + —4, + —4 - — + - + — <(4), 

2304 1152 1024 768 4096 576 34560 2 

1 1 1 Q ^2 1 q^4 „4 



19 4 4 


13 4 


+ 15 Is- 


768 


21 


25tt 2 1 


2304 


1152 


1024 3 


4096 


+ 576 


- 1 4 


+ 1 I 


13 


^ 2 


13tt 4 


+ y*?(4) 


1152 


1152 


2 55296 


~ 452 ~ 


30240 


57ef + 


7 44 


- 1 is- 


3^ 4 + 
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17tt 2 


4331tt 4 
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-$ 2 '(4), 

4d 61 - 37 - 105- 59 - 17tt 2 4331^ 4 , , , 

7T 4 B = -TTTTT^i - ^7^2 - —4 ~ ^4 - 4 og/|n + tt 4 $^(4) - 4tt 4 $^(4), 

K *B, = —4 - —4 - — + — + 33l7r4 + 7r 4 $''(4) - 6^ 4 $'"(4) 

1 96 1 96 2 4608 + 288 + 30240 + ll ' 1 1 h 

tt 4 ^ = —I - —4 —4 + —4 - — + 40971 " 4 - tt 4 $' f4) - 47T 4 $"(4) 

2 576 1 288 2 1536 3 + 384 4 48 + 4320 21 J 21 j ' 

, 23 - 11 - 45 - 37- 13tt 2 421tt 4 4 , , d s , .. , s 

* = 288* + 144^ 2 + 128 3 + 96 + ~L8 720" + ^ ~ ^ ^ + 
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144 18 384 96 18 720 



144 18 384 J 96 18 720 w 2V ' 2W ' 

4n 1 - 1 - 157 52tt 2 253tt 4 
io 



18 18 1728 27 840 ' 

In order to compare the two approaches we have not included the dependence on l/F® in the parameters In 
this way the quantities bi that we have used in Eqs. (0) stand for: 
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-16 fcl -^6- fc2+r5 ' 
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K1 + r R , 

48 1 Qfi 
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" 1 /j 



1927T 4 fj, 

and /_t is the renormalization mass scale. 
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